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Executive  Summary 

All  military  organizations  depend  on  the  reliable  performance  of  repairable  systems  for 
the  successful  completion  of  missions.  Due  to  limitations  in  maintenance  resources,  a 
maintenance  manager  must  decide  how  to  allocate  available  resources.  This  allocation  falls 
within  the  domain  of  selective  maintenance.  Selective  maintenance  is  defined  as  the  process  of 
identifying  the  subset  of  maintenance  activities  to  perform  from  a  set  of  desired  maintenance 
actions.  Previously,  researchers  have  developed  a  class  of  mathematical  models  that  can  be  used 
to  identify  selective  maintenance  decisions  for  the  following  scenario  -  A  system  has  just 
completed  a  mission  and  will  begin  its  next  mission  soon.  Maintenance  cannot  be  performed 
during  missions;  therefore,  the  decision-maker  must  decide  which  components  to  maintain  prior 
to  the  next  mission.  The  selective  maintenance  models  considered  to  date  treat  decision-making 
for  binary-state  systems,  i.e.  all  components,  subsystems,  and  the  system  itself  are  assumed  to  be 
either  functioning  or  failed  at  any  point  in  time. 

The  primary  objective  of  this  project  is  to  develop  a  modeling-based  methodology  for 
managing  selective  maintenance  decisions  when  multiple  (more  than  two)  system  states  are 
possible.  First,  the  research  literature  for  selective  maintenance  and  multi-state  analysis  is 
presented.  Then,  we  define  a  scenario  in  which  systems  in  various  states  of  maintenance  need 
must  perform  a  number  of  different  missions.  For  this  scenario,  we  formulate  a  non-linear 
mathematical  program,  and  we  explore  three  solution  procedures  for  the  optimization  problem. 
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1.  Introduction 


All  military  organizations  depend  on  the  reliable  performance  of  repairable  systems  for 
the  successful  completion  of  missions.  The  use  of  mathematical  modeling  for  the  purpose  of 
modeling  repairable  systems  and  designing  optimal  maintenance  policies  for  these  systems  has 
received  an  extensive  amount  of  attention  in  the  literature.  Unfortunately,  the  vast  majority  of 
this  work  ignores  potential  limitations  on  the  resources  required  to  perform  maintenance  actions. 
This  shortcoming  has  motivated  the  development  of  models  for  selective  maintenance,  the 
process  of  identifying  the  subset  of  actions  to  perform  from  a  set  of  desirable  maintenance 
actions.  Previously,  we  have  developed  a  class  of  mathematical  models  that  can  be  used  to 
identify  selective  maintenance  decisions  for  the  following  scenario  -  A  system  has  just 
completed  a  mission  and  will  begin  its  next  mission  soon.  Maintenance  cannot  be  performed 
during  missions;  therefore  the  decision-maker  must  decide  which  components  to  maintain  prior 
to  the  next  mission.  The  selective  maintenance  models  formulated  to  date  are  based  on  the 
assumption  of  binary  (functioning,  failed)  component,  subsystem  and  system  status.  As  a  result, 
mission  reliability  is  used  as  the  objective  function  in  the  resulting  mathematical  programming 
models.  We  wish  to  improve  upon  this  approach  in  two  ways.  First,  it  may  be  more  realistic  to 
classify  component  status  using  more  than  two  discrete  levels  (if  not  some  continuous  measure). 
This  implies  multi-state  measures  of  subsystem  and  system  status  as  well.  Second,  the 
performance  of  a  military  system  typically  can  be  measured  using  several  measures  in  addition  to 
mission  reliability.  All  these  performance  measures  are  functions  of  the  status  of  the 
components.  The  primary  objective  of  this  project  is  to  develop  multi-state  selective  maintenance 
models  that  incorporate  multi-state  component  status  and  multiple  measures  of  system 
performance. 
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The  activities  required  to  achieve  the  objective  of  this  project  are  applied  to  a  set  of 
systems  utilized  by  the  US  Air  Force.  First,  we  define  the  system  structure  and  appropriate  status 
measures  for  each  component  in  the  system.  Second,  we  identify  the  resources  consumed  by 
maintenance  actions,  the  impact  on  component  status  of  each  potential  maintenance  action,  and 
the  quantity  of  each  resource  consumed  by  each  maintenance  action.  Third,  we  identify  the 
relevant  measures  of  mission  performance  and  develop  functions  which  capture  these  measures 
in  terms  of  the  component  status  values.  Fourth,  we  develop  a  mathematical  formulation  of  the 
selective  maintenance  problem.  Finally,  we  develop  solution  procedures  for  solving  the  selective 
maintenance  problem.  We  define  enumerative  solution  strategies  for  smaller  problems  and 
heuristic  strategies  for  larger  problems. 
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2.  Research  Literature  Review 


In  this  section,  a  review  of  the  relevant  research  literature  is  presented.  We  begin  with  a 
summary  of  the  selective  maintenance  literature.  Then,  we  conclude  with  a  summary  of  the 
literature  related  to  multi-state  reliability  analysis. 

2.1  Selective  Maintenance 

This  project  builds  upon  the  body  of  knowledge  in  selective  maintenance.  Selective 
maintenance  falls  within  the  domain  of  maintenance  modeling  and  optimization.  The  use  of 
mathematical  modeling  for  the  purpose  of  modeling  repairable  systems  and  designing  optimal 
maintenance  policies  for  these  systems  has  received  an  extensive  amount  of  attention  in  the 
literature  [5, 6, 7,  8, 10, 12, 13]. 

The  original  study  in  selective  maintenance  was  performed  by  Rice  et  al.  [9].  They  define 
a  system  that  must  complete  a  series  of  missions  where  maintenance  is  performed  only  during 
finite  breaks  between  missions.  Due  to  the  limited  maintenance  time,  it  may  not  be  possible  to 
repair  all  failed  components  before  the  next  mission.  A  nonlinear,  discrete  selective  maintenance 
optimization  model  is  developed  which  is  designed  to  maximize  system  reliability  for  the  next 
mission.  The  numbers  of  components  to  repair  are  the  decision  variables,  and  the  limitation  on 
maintenance  time  serves  as  the  primary  functional  constraint.  Due  to  the  complexity  of  the 
model,  total  enumeration  is  the  recommended  solution  procedure.  Given  that  total  enumeration  is 
ineffective  for  large  scenarios,  a  heuristic  selective  maintenance  procedure  is  developed. 

Cassady  et  al.  [1,  2]  extend  the  work  of  Rice  et  al.  [9]  in  several  ways.  First,  more 
complex  systems  are  analyzed.  Specifically,  systems  are  comprised  of  independent  subsystems 
connected  in  series  with  the  individual  components  in  each  subsystem  connected  in  any  fashion. 
Next,  the  selective  maintenance  model  is  extended  to  consider  the  case  where  both  time  and  cost 
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are  constrained.  This  leads  to  the  development  of  three  different  selective  maintenance  models. 
These  models  include  maximizing  system  reliability  subject  to  both  time  and  cost  constraints; 
minimizing  system  repair  costs  subject  to  a  time  constraint  and  a  minimum  required  reliability 
level;  and  minimizing  total  repair  time  subject  to  both  cost  and  reliability  constraints. 

Cassady  et  al.  [3]  extend  the  work  of  Rice  et  al.  [9]  in  two  other  ways.  First,  system 
components  are  assumed  to  have  Weibull  life  distributions.  This  assumption  permits  systems  to 
experience  an  increasing  failure  rate  (IFR)  and  requires  monitoring  of  the  age  of  components. 
Second,  the  selective  maintenance  model  is  formulated  to  include  three  maintenance  actions: 
minimal  repair  of  failed  components,  replacement  of  failed  components,  and  preventive 
maintenance. 

Chen  et  al.  [4]  extend  the  work  of  Rice  et  al.  [9]  and  Cassady  et  al.  [1]  by  considering 
systems  in  which  each  component  and  the  system  may  be  in  K  +  1  possible  states,  0,  1,  ...  ,  K. 
They  use  an  optimization  model  to  minimize  the  total  cost  of  maintenance  activities  subject  to  a 
minimum  required  system  reliability. 

Schneider  and  Cassady  [10]  formulate  an  optimization  model  to  extend  the  work  of  Rice 
et  al.  [9]  by  defining  a  selective  maintenance  model  for  a  set  of  systems  that  must  perform  a  set 
of  missions  with  system  maintenance  performed  only  between  sets  of  missions.  Three  models  are 
formulated.  The  first  model  maximizes  the  probability  that  all  systems  within  the  set  successfully 
complete  die  next  mission,  where  as  the  second  model  minimizes  the  variable  cost  associated 
with  maintenance.  A  special  case  of  the  second  model  allows  the  user  to  maximize  the  expected 
value  of  the  number  of  successful  missions  in  the  next  set.  The  third  model  permits  cancellation 
of  a  mission  based  on  costs  associated  with  the  risk  of  failure. 
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2.2  Multi-State  Analysis 

Recently  a  great  deal  of  attention  has  been  given  to  assessing  and  optimizing  multi-state 
systems.  For  a  comprehensive  literature  review  of  multi-state  systems,  see  Lisnianski  and  Levitin 
[14].  However,  all  of  this  research  focuses  on  systems  with  coherent  states.  States  exhibit 
coherency  when  there  is  a  spectrum  of  states  with  perfectly  functioning  on  one  end  and 
completely  failed  at  the  other  and  each  state  in  between  is  a  uniform,  incremental  gradation  or 
degradation  towards  being  perfectly  functioning  or  completely  failed.  An  example  would  be  tire 
tread  wear.  The  two  ends  of  the  tread  wear  spectrum  are  brand  new  and  completely  bald  tires; 
these  states  may  be  denoted  1 0  and  0  respectively.  The  states  in-between  may  be  represented  by 
the  whole  numbers  {1, 2, ... ,  9}  where  each  state  corresponds  to  the  percentage  of  tread 
remaining:  state  9:  90%,  state  8:  80%,  etc.  Coherent  states  imply  that  10  is  better  than  9,  9  is 
better  than  8,  etc.,  and  that  the  difference  in  states  is  uniform  whereby  the  difference  between 
states  8  and  9  is  the  same  as  the  difference  between  states  2  and  3. 

This  research  deals  with  non-coherent  states  that  are  descriptive  or  qualitative  in  nature. 
The  difference  in  the  states  are  not  necessarily  incrementally  increasing  or  decreasing,  but  denote 
which  missions  the  system  is  capable  of  performing  rather  than  the  reliability  of  the  system  or 
some  other  performance  metric.  The  state  of  the  system  is  denoted  by  a  vector  of  the  binary 
status  of  all  of  the  subsystems  of  the  systems.  Thus  far  no  research  has  been  conducted  in 
optimizing  systems  with  multiple,  non-coherent  states. 
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3.  Problem  Statement 


3.1  Scenario  of  Interest 

Consider  the  following  scenario.  There  are  q  systems,  each  comprised  of  m  independent 
subsystems  that  are  idle  and  available  for  maintenance.  The  state  of  system  /  is  denoted  aj  =  (an, 
an,  ■  ■  ■ ,  aim)  where  ay  denotes  the  amount  of  time  required  to  bring  subsystem  j  of  system  i  into 
a  properly  operating  condition  where  ay  e  {0,1 , ,<pj).  Table  1  displays  an  example  a  vector 
(q  =  24,  m  =  41)  that  will  be  used  as  data  for  the  numerical  example  we  will  be  presenting 
throughout  this  paper. 

Some  maintenance  actions  require  spare  parts  or  other  resources  that  are  not  readily 
available.  We  capture  this  by  assigning  a  ready  time  to  each  subsystem.  The  ready  time  of 
subsystem  j  in  system  i,  denoted  py,  is  the  time  at  which  these  resources  are  available  and 
maintenance  on  the  subsystem  can  begin.  The  ready  time  of  system  i  is  denoted  by  the  vector  pj 
=  (pn,  pa, . . . ,  Pim).  Table  2  presents  an  example  p  vector  for  our  numerical  example. 

3.2  Mission  Profile 

There  are  n  future  missions  planned  where  n  <  q.  Mission  k  requires  some  subset  of  the 
systems  to  be  operational.  Therefore,  the  state  of  mission  k  is  denoted  =  (ski,  Sk2,  ■  •  .  ,  Shn) 
where 


fl  if  mission  k  requires  subsystem  j 


kj  1 0  otherwise 


An  example  s  vector  (k  =  24)  for  our  numerical  example  is  presented  in  Table  3. 
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Table  1 .  Example  a  Vector 


System 


Subsystem 

1 

2 

3 

4 

5  6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

Air  Frame 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Crew  Station  System 

0 

0 

6 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

2 

0 

Landing  Gear  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Flight  Control  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Turbofan  Power  Plant  (PW  Engines) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

Aux  Power  Plant/JFS 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

7 

0 

0 

0 

0 

0 

0 

Turbofan  Power  Plant  (GE  110) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Environmental  Control  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 

50 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Electrical  Power  Supply 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Exterior  Lighting 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

Interior  Lighting 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Hydraulic/Pneumatic  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Fuel  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

Oxygen  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Fire  Detection  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

Overheat  DetectionSystem 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Flight  Instrumentals 

Malfunction  Analysis  &  Recording 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

1 

Equipment  (CFSDR) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

HF  Communications 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

VHF  Communications 

0 

0 

0 

0 

0 

0 

48 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

UHF  Communications 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

Interphone  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

IFF 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

9 

0 

92 

0 

0 

Improved  Data  Modem  (IDM)/Situation 

Awareness  data  link  (SADL) 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Radio  Navigation 

0 

0 

2 

p 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Global  Positioning  System  (GPS) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Fire  Control  System 

0 

2 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

Airborne  Video  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

0 

20 

0 

0 

0 

0 

0 

0 

0 

Data  Transfer  Unit 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Radar  Altimeter  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

75 

0 

0 

0 

0 

Targeting  Pod  (GTP)  System 

0 

29 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Navigation  Pod  (VP)  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

2 

Harm  Targeting  System  (HTS) 

0 

21 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.0 

Weapons  Delivery  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Gun  System 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

Electronic  Counter  Measures 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Radar  Warning  Receiver 

0 

6 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Chaff/Flare  Disp  System 

8 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Emergency  Equipment 

1 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Smoke  Generation  System 

0 

0 

0 

0 

0 

0 

31 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Explosive  Devices  and  Components 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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Table  2.  Example  q_  Vector 


System 


Subystem 

1 

2  3 

4 

5 

6 

7 

8  9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

Air  Frame 

0 

0  0 

0 

0 

0 

96 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

Crew  Station  System 

2 

0  0 

0 

0 

0 

96 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Landing  Gear  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

4 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Flight  Control  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

0 

0 

0 

0 

0 

0 

Turbofan  Power  Plant  (PW  Engines) 

0 

0  0 

% 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

Aux  Power  Plant/JFS 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

96 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Turbofan  Power  Plant  (GE  110) 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Environmental  Control  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Electrical  Power  Supply 

0 

0  0 

0 

0 

0 

0 

4  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

Exterior  Lighting 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Interior  Lighting 

0 

0  0 

0 

0 

96 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

Hydraulic/Pneumatic  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Fuel  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Oxygen  System 

0 

0  0 

96 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Fire  Detection  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Overheat  DetectionSystem 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

4 

4 

0 

Flight  Instrumentals 

Malfunction  Analysis  &  Recording 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

0 

0 

0 

0 

0 

0 

Equipment  (CFSDR) 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

HF  Communications 

0 

0  0 

0 

0 

96 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

VHF  Communications 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

0 

0 

UHF  Communications 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Interphone  System 

0 

0  0 

0 

0 

0 

0 

0  0 

4 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

IFF 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Improved  Data  Modem  (IDM)/Situation 

Awareness  data  link  (SADL) 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

Radio  Navigation 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Global  Positioning  System  (GPS) 

0 

0  0 

0 

0 

0 

0 

0  0 

4 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Fire  Control  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Airborne  Video  System 

2 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Data  Transfer  Unit 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

Radar  Altimeter  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

0 

Targeting  Pod  (GTP)  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Navigation  Pod  (VP)  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Harm  Targeting  System  (FITS) 

0 

0  4 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

Weapons  Delivery  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

96 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Gun  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

Electronic  Counter  Measures 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

% 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Radar  Warning  Receiver 

0 

0  0 

0 

96 

96 

0 

0  0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

Chaff/Flare  Disp  System 

0 

0  0 

0 

0 

0 

0 

96  0 

4 

96 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

Emergency  Equipment 

0 

0  0 

0 

0 

0 

0 

0  0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Smoke  Generation  System 

0 

0  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

96 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Explosive  Devices  and  Components 

0 

0  0 

0 

0 

0 

0 

2  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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Table  3.  Example  s  Vector 


Subsystem 
Air  Frame 

Crew  Station  System 

Landing  Gear  System 

Flight  Control  System 

Turbofan  Power  Plant  (PW  Engines) 

Aux  Power  Plant/JFS 

Turbofan  Power  Plant  (GE  1 10) 

Environmental  Control  System 

Electrical  Power  Supply 

Exterior  Lighting 

Interior  Lighting 

Hydraulic/Pneumatic  System 

Fuel  System 

Oxygen  System 

Fire  Detection  System 

Overheat  DetectionSystem 

Flight  Instrumentals 

Malfunction  Analysis  &  Recording 

Equipment  (CFSDR) 

HF  Communications 
VHF  Communications 
UHF  Communications 
Interphone  System 
IFF 

Improved  Data  Modem  (IDM)/Situation 
Awareness  data  link  (SADL) 

Radio  Navigation 

Global  Positioning  System  (GPS) 

Fire  Control  System 
Airborne  Video  System 
Data  Transfer  Unit 
Radar  Altimeter  System 
Targeting  Pod  (GTP)  System 
Navigation  Pod  (VP)  System 
Harm  Targeting  System  (HTS) 

Weapons  Delivery  System 
Gun  System 

Electronic  Counter  Measures 
Radar  Warning  Receiver 
Chaff/Flare  Disp  System 
Emergency  Equipment 
Smoke  Generation  System 
Explosive  Devices  and  Components 


Mission 


FSL  ADC  ASC  ASY 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
till 
1111 
1111 
1111 
1111 
1111 

10  0  0 
110  0 
1111 
1111 
1111 
1111 

1111 
1111 
1111 
1111 
10  0  0 

1111 
10  10 

10  10 

10  10 

10  0  0 

11  1  1 

1111 
1111 
1111 
1111 
1111 
10  0  0 

1111 


ASN  DSP  TNG  TST 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
1111 
11  1  1 

0  0  0  0 

0  0  11 

1111 
1111 
11  1  1 

1111 

1111 
1111 
1111 
1111 
0  0  0  0 

1111 
10  11 
10  11 
10  11 
0  10  0 

1111 
0  110 
110  1 
1111 
1111 
1111 
0  0  0  0 

1111 
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3.3  Decision  Making  I 

It  must  be  determined  which  system  should  be  assigned  to  each  mission.  The  decision 
variable 


{1  if  system  i  is  assigned  to  mission  k 
0  otherwise 

is  used  to  denote  the  assignments.  It  will  be  assumed  that  every  mission  is  assigned  a  system  and 
no  system  is  assigned  more  than  one  mission.  Once  the  assignments  have  been  made,  the  total 
time  required  for  maintenance  related  to  mission  k  is 


=2Xfe*  • 9. ,) 

<=i 

This  is  generally  referred  to  as  the  processing  time  in  the  scheduling  literature.  It  is  the  sum  over 
all  systems  of  the  dot  product  of  the  state  vector  and  the  condition  vector,  given  that  the 
assignment  is  made.  The  condition  of  a  subsystem  may  be  negated  if  the  mission  does  not 
require  the  subsystem.  Likewise  the  ready  time  for  maintenance  related  to  mission  k  is 


max 


For  this  problem  we  assume  that  the  waiting  time  for  parts  or  various  delays  occurs  concurrently. 
Therefore,  the  ready  time  of  the  mission  is  just  the  highest  ready  time  of  subsystem  required  by 
the  mission. 

3.4  Decision  Making  II 

Once  the  assignments  of  systems  to  missions  are  made,  maintenance  crews  must  perform 
the  maintenance  ( y  denotes  the  number  of  crews  available  at  one  time).  We  assume  that  a  crew 
(1)  works  on  one  system  at  a  time,  (2)  works  on  a  system  after  it  is  “ready”,  and  (3)  works  on  a 
system  until  all  maintenance  is  finished  (i.e.  no  preemptions). 
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For  each  mission,  a  decision  must  be  made  as  to  when  maintenance  will  begin  and  by 


which  crew  the  maintenance  will  be  performed.  A  second  decision  variable 


if  mx  for  mission  k  is  initiated  by  crew  /  at  time  t 
otherwise 


denotes  if  maintenance  for  mission  k  will  be  initiated  by  crew  /  at  time  t.  Time  is  considered  to 
be  a  discrete  index  where  /  e  {0, 1, ... ,  Cmax  -l}  .  The  time  span  (Cmax)  is  defined  as 

max(Ci,....,C„)  where  Cj  is  the  completion  time  of  job  j.  This  is  equivalent  to  the  completion  time 
of  the  final  job  in  the  system  [15].  Since,  in  our  scenario,  Cmax  is  a  function  of  decision  variables, 
we  will  place  the  following  upper  bound  on  Cmax 

Cm«x=XZk+A/) 

i'=l  7=1 

The  following  assumptions  will  be  made  in  the  second  decision  making  process:  (1)  every 
mission  gets  a  crew,  (2)  maintenance  cannot  be  started  before  the  mission  is  “ready”,  and  (3)  a 
crew  cannot  work  on  two  systems  at  the  same  time. 

3.5  Objective 

The  overall  objective  of  this  problem  is  to  minimize  the  total  weighted  completion  time 
(TWCT)  of  maintenance  on  all  missions.  Weighted  completion  time  is  also  referred  to  as 
weighted  flow  time  in  the  literature.  The  problem  is  actually  a  combination  of  an  assignment 
problem  and  a  scheduling  problem.  Systems  must  be  assigned  to  missions  such  that  the  total 
maintenance  time  and  delay  time  due  to  ready  times  of  the  assignment  corresponds 
simultaneously  with  minimizing  the  TWCT  of  the  scheduled  maintenance. 

The  scheduling  problem  is  denoted  Pm  |  ry  |  ^  wyC  •  using  the  notation  of  Lawler  et  al. 
[16].  The  scheduling  problem  is  considered  to  have  identical  parallel  machines  {Pm)  with  ready 
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times  (rj)  and  a  TWCT  objective.  The  crews  performing  maintenance  are  considered  to  be  the 
identical  parallel  machines  as  it  will  be  assumed  that  all  maintenance  crews  contain  the  same 
number  of  people  and  can  perform  maintenance  at  the  same  pace.  We  are  assigning  the  various 
missions  different  weights  (wj)  or  importance  because  it  may  be  desired  for  more  important 
missions  to  have  an  earlier  completion  time  (Q).  A  larger  weight  implies  more  importance. 
Weights  for  our  numerical  example  are  presented  in  Table  4.  The  Pm  |  ry  |  ^w,C,  is  strongly 
NP-hard. 


Table  4.  Example  Weights  (w*) 


Mission  Type 

FSL 

8 

ADC 

6 

ASC 

10 

ASY 

5 

ASN 

9 

DSP 

3 

TNG 

1 

TST 

2 

12 


4.  Motivation 


This  research  is  motivated  by  a  real-world  Air  Force  scenario  in  which  planes  can  have 
multiple  qualitative  states.  The  system  under  consideration  is  the  F-16  A/B/C/D  due  to  our 
experience  with  them  at  Hill  AFB,  UT.  There  are  41  subsystems  for  the  F-16  as  detailed  by  the 
mission  essential  subsystem  list  (MESL)  (AFI121-103  ACCSUP1).  There  are  also  eight  possible 
missions  for  an  F-16: 

1 .  FSL  -  Full  System  List 

2.  ADC  -  Air  Defense,  Conventional 

3.  ASC  -  Air  to  Surface,  Conventional 

4.  ASY  -  Air  Superiority 

5.  ASN  -  Air  to  Surface,  Nuclear 

6.  DSP  -  Defense  Suppression 

7.  TNG  -  Training 

8.  TST  -  Testing 

Mission  FSL  denotes  that  all  subsystems  are  required.  The  remaining  seven  missions  require 
some  combination  of  the  41  subsystems,  but  not  all.  The  qualitative  state  of  each  plane  is  which 
missions  it  is  capable  of  performing. 


13 


5.  Mathematical  Formulation 


Sets  and  Indexes 

i  =  system  number  (/  =  1,2, ...  ,q) 
j  =  subsystem  number  (j  —  1,2, 
k  =  mission  number  (k  =  1,2, . . .,  n) 
l  s  crew  number  (/  =1,2,...,  y) 

/  =  time  (t  ~  0,  1 , ... ,  Cmax  ~1 ) 

Parameters 

ay  =  amount  of  time  required  to  bring  subsystem  j  of  system  i  into  a  properly  operating  condition 

{1  if  mission  k  requires  subsytem  j 

0  otherwise 

py  =  ready  time  of  subsystem  j  in  system  i 
M>k  =  importance  (weight)  of  mission  k 


Decision  Variables 


x 


ik 


1  if  system  /  is  assigned  to  mission  k 
0  otherwise 


yu<  = 


1  if  maintenance  for  mission  k  is  initiated  by  crew  /  at  time  / 
0  otherwise 
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5.1  Mathematical  Program 


minimize 


subject  to: 


£wki 

*=i 


ft'  • 


;= 1  t=o 


1=1  1  jm\ 


yu, 


lx =1 

P*sl 

/=1  1=0 
n  /-I 

E  E^--1 

*=1  /'=max(»-pj,0) 

c*  ~Pk^rk 

9 


^=Eu&}^h 


k  =  1,2 . « 

/  =  1,2,. ..,9 

A:  =  1,2,...,  rt 

/  =  1,2 . y 

t  =  0,1, . . . ,  Cmax  —1 
A;  =  1,2, _ , « 

k  =  l,2,...,n 


/= 1  /=o 


9 

(  m  ) 

*  +  Z 

Xik  E  a0  Skj 

/=1 

\  7=1  ) 

_ 

«  I  7=1  J 

cm»=i!Ek+^) 

i=l  7=1 

*/*=0orl 


yklt  =  Oor  1 


k  =  1,2,...,» 


j  =  1,2,...,<7 

k  =  1,2, ...,w 
£  =  1,2,...,« 

/  =  l,2,...,y 
t  =  0,l,...,Cmax  -1 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 


The  above  formulation  represents  the  non-linear  binary  program  for  the  multi-state 
selective  maintenance  problem.  The  formulation  is  constructed  using  time-indexed  variables. 
The  objective  function  minimizes  TWCT  where  the  completion  time  of  mission  k  is  given  by  (7) 
in  the  list  of  constraints.  Constraint  (1)  says  that  every  mission  is  assigned  a  system,  and  likewise 
(2)  ensures  that  no  system  gets  more  than  one  mission.  Constraint  (3)  ensures  that  every  mission 
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is  assigned  a  crew.  Constraint  (4)  does  not  allow  preemption  or  the  ability  for  a  crew  to  work  on 
two  systems  at  the  same  time.  Constraint  (5)  enforces  the  ready  times  of  the  missions  by 
ensuring  that  the  completion  time  of  the  mission  minus  its  processing  time  is  greater  than  or 
equal  to  the  ready  time  of  the  mission.  Constraints  (6-7)  set  the  ready  time  and  completion  time 
of  the  missions,  respectively  as  previously  defined  in  the  problem  statement.  Constraints  (10-11) 
enforce  the  binary  restriction  of  the  decision  variables. 

As  previously,  mentioned  the  scheduling  problem  alone  is  strongly  NP-hard  and  therefore 
the  entire  problem  is  strongly  NP-hard.  With  the  combination  assignment  problem  and 
scheduling  problem  the  formulation  is  non-linear  both  in  the  objective  function  and  in  the 
constraints.  Non-linearity  occurs  in  the  objective  function  when  the  two  decision  variables  are 
multiplied  together  when  the  summations  are  expanded.  This  specific  case  of  non-linearity  is 
known  as  a  quadratic  objective  function.  Non-linearity  is  also  introduced  in  constraint  (4)  when 
the  lower  bound  of  the  inner  summation  contains  the  variable  pk  which  is  a  function  of  the 
decision  variable  *,*. 
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6.  Solution  Procedure 


Three  solution  approaches  for  the  multi-state  selective  maintenance  problem  were 
investigated.  First,  the  problem  was  solved  using  a  total  enumeration  procedure.  Second,  a 
combination  of  a  heuristic  and  commercial  solver  was  utilized.  Finally,  a  dispatching  rule  was 
tested. 

6.1  Total  Enumeration 

The  total  enumeration  procedure  was  coded  using  Visual  Basic  for  Applications  (VBA) 
running  behind  a  Microsoft  Excel  interface.  It  is  an  exponential  time  procedure  that  first 
enumerates  all  of  the  possible  assignments  and  for  each  assignment  enumerates  all  possible 
schedules  and  reports  the  minimum  objective  value  found. 

Due  to  the  computational  complexity  of  the  multi-state  selective  maintenance  problem, 
only  extremely  small  problem  instances  can  be  solved  in  a  reasonable  amount  of  time  using  total 
enumeration.  The  number  of  iterations  required  to  enumerate  the  multi-state  selective 
maintenance  problem  is  1 02'7.  Thus,  for  a  three  mission,  three,  system,  two  crews  scenario,  one- 
trillion  iterations  are  required.  A  Pentium  IV  2.0  GHz  with  2.0  GB  of  RAM  took  weeks  to 
enumerate  this  small  problem  instance.  Obviously,  a  scenario  this  small  is  not  useful  for 
application  by  the  Air  Force  and  even  if  it  was,  this  is  not  a  practical  computation  time  for 
implementation.  Therefore,  other  solutions  approaches  must  be  explored. 

6.2  Heuristic  and  Commercial  Solver 

A  commercial  solver  has  been  utilized  to  facilitate  the  solution  process.  Unfortunately, 
the  formulation  is  a  non-linear  program  which  makes  it  incapable  of  being  solved  by  most  almost 
all  commercial  solver  packages.  Some  packages  are  capable  of  solving  quadratic  programs  as  our 
problem  formulation  contains  a  quadratic  objective  function.  However,  the  non-linear  constraint 
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(4)  eliminates  the  quadratic  properties  of  the  problem.  Therefore,  we  have  implemented  a 
heuristic  in  combination  with  a  commercial  solver  to  obtain  solutions  to  the  multi-state  selective 
maintenance  problem.  The  heuristic  simply  breaks  the  multi-state  selective  maintenance  problem 
up  into  two  separate  subproblems:  an  assignment  problem  and  a  scheduling  problem.  The 
assignment  problem  will  be  solved  to  optimality  using  total  enumeration  for  a  ready-time-based 
weighted  shortest  processing  time  (WSPT)  first  objective  function.  The  heuristic  is  based  on  the 
famous  rule  in  scheduling  theory  that  jobs  that  are  ordered  in  decreasing  order  of  Wj  /  pj  is 
optimal  for  the  single  machine  TWCT  (1 1|  ^w.C, )  problem.  Our  problem  deviates  from  this 

simplified  problem  in  that  it  is  a  more  complex  version  containing  ready-times  and  multiple 
parallel  machines.  However,  the  modified  objective  function  will  still  serve  to  reduce  the  TWCT 
objective.  The  heuristic  formulates  assignments  of  systems  to  missions  by  making  the 
assignment  with  the  highest 

w*  _ w* _ 

Pt+^k  !>*& 

And  continuing  in  this  fashion  until  all  missions  have  been  assigned  a  system  and  not  more  than 
one  system  is  assigned  a  mission.  The  heuristic  is  denoted  a  read-time  based  WSPT  because  of 
the  addition  of  rk  in  the  denominator.  It  is  desired  to  have  the  missions  with  the  highest  wk  / pk  be 
scheduled  as  soon  as  possible  so  the  ready  time  is  placed  in  the  denominator  to  have  it 
minimized.  The  square  root  of  the  ready-time  is  taken  so  it  does  not  dominate  the  processing 
time.  Ready-times  will  often  be  several  orders  of  magnitudes  higher  than  the  processing  time  as 
we  are  measuring  time  in  hours  and  ready  times  may  be  several  days.  The  heuristic  for  the 
assignment  problem  was  coded  in  VBA  and  runs  in  a  fraction  of  a  second.  Table  6  demonstrates 
the  assignment  made  by  the  heuristic  for  our  numerical  example  data  we  have  presented 
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previously.  For  our  example,  consider  the  distribution  of  the  number  of  each  missions  type 
presented  in  Table  5.  The  assignment  column  in  Table  6  shows  which  system  is  assigned  to 
which  mission  type.  The  job  formulation  column  displays  the  resulting  job  parameters  of  the 
assignment.  After  the  assignment  has  been  made,  each  job  will  have  an  associated  weight  (wj), 
ready  time  (rj),  and  processing  time  (pj) 

Table  5.  Numerical  Example  Mission  Type  Distribution 


Mission  Type 

# 

FSL 

7 

ADC 

3 

ASC 

4 

ASY 

5 

ASN 

1 

DSP 

2 

TNG 

2 

TST 

0 

19 


Table  6.  Numerical  Example  Heuristic  Assignment 


Assignment 

|  Job  Formulation 

System 

Mission 

# 

wi 

n 

Pi 

13 

FSL 

1 

1 

96 

0 

21 

ADC 

2 

1 

96 

0 

19 

FSL 

3 

10 

0 

1 

10 

FSL 

4 

9 

4 

1 

9 

FSL 

5 

10 

0 

4 

23 

FSL 

6 

6 

4 

2 

17 

ASY 

7 

10 

2 

7 

12 

FSL 

8 

10 

4 

8 

1 

ASC 

9 

8 

2 

9 

3 

ASN 

10 

8 

4 

11 

4 

FSL 

11 

8 

96 

1 

5 

ASC 

12 

8 

96 

3 

11 

TNG 

13 

8 

96 

4 

8 

ASY 

14 

8 

96 

2 

6 

ASY 

15 

6 

96 

4 

2 

ADC 

16 

6 

0 

8 

16 

ASC 

17 

3 

0 

11 

24 

ASY 

18 

3 

96 

1 

22 

ASC 

19 

5 

96 

92 

20 

ASY 

20 

5 

2 

21 

18 

TNG 

21 

5 

96 

15 

15 

DSP 

22 

5 

4 

53 

14 

ADC 

23 

5 

96 

10 

7 

DSP 

24 

5 

96 

51 

The  heuristic  solution  to  the  assignment  problem  will  then  be  used  as  an  input  to  the 
scheduling  problem  which  is  formulated  in  AMPL  and  solved  by  ILOG’s  CPLEX  9.0.  CPLEX 
can  obtain  the  optimal  solution  to  the  scheduling  problem,  however,  a  global  optimal  solution  for 
the  aggregated  problem  cannot  be  guaranteed,  because  the  solution  of  the  assignment  problem  is 
not  necessarily  the  assignment  for  the  global  optimum  of  the  entire  problem.  The  mathematical 
formulation  of  the  scheduling  problem  is  a  subset  of  the  entire  multi-state  selective  maintenance 
formulation  and  is  presented  below.  The  AMPL  formulation  is  contained  in  Appendix  A. 
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Sets  and  Indexes 

k  =  mission  number  (k  =  1, 2, . . . ,  n) 
l  =  machine  number  (l  =  1,2, ...  ,y) 
t  =  time  (/  =  0, 1, . . . ,  Cmax  -1) 


Parameters 

Pk  =  processing  time  of  job  A: 

Wk  =  weight  of  job  k 

Decision  Variable 

{1  if  job  k  starts  processing  at  time  t  by  crew  1 
0  otherwise 


minimize 


*=1  /=1  1=0 


subject  to: 


£  c2>»  = 


1 


/=1  /= 0 
n  /-I 


Z  Z***'  -1 

*=1  t'=max\t-pj  o ) 

££■ 

/= 1  r=0 


*U,  =  Sk 


Xyi  =  0  or  1 


£  =  1,2,...,/! 

1  =  1,2,...,;' 
t  =  0,l,...,Cmax-l 

k  =  1,2,...,« 

A:  =  l,2,...,n 
A:  =  l,2,...,w 

/  =  0,1, . . . ,  Cmax  —  1 


The  optimal  crew  assignments  and  schedule  obtained  by  CPLEX  for  our  numerical 
example  is  displayed  in  Table  7.  The  TWCT  objective  function  for  our  numerical  example  given 
by  the  heuristic/CPLEX  solution  method  is  7,992. 
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Table  7.  Numerical  Example:  CPLEX  Crew  Assignments  and  Schedule 


Crew  # 

Job# 

Start  Time  (t) 

1 

9 

2 

12 

96 

19 

99 

2 

7 

2 

20 

9 

15 

96 

3 

16 

0 

22 

8 

14 

96 

21 

98 

4 

3 

0 

4 

4 

8 

5 

18 

96 

24 

97 

5 

5 

0 

6 

4 

10 

6 

11 

96 

23 

97 

6 

17 

0 

13 

96 

1 

n/a 

2 

n/a 

6.3  Dispatching  Rule 

The  dispatching  rule  is  a  simplified  version  of  the  heuristic  presented  above  without  the 
CPLEX  optimization.  The  dispatching  rule  is  designed  to  be  easily  calculated  by  hand  and 
implemented  in  the  field  without  the  aid  of  a  computer.  Missions  are  simply  “dispatched”  or 
scheduled  in  order  of  the  decreasing  objective  function 


m 

I X+-, 

H 

>=! 
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For  each  system  i  the  total  processing  times  of  all  m  subsystems  are  summed  and  added  to  the 
square  root  of  the  sum  of  all  of  the  ready  times  of  each  subsystem  j  on  system  /.  Missions  are 
dispatched  by  taking  the  mission  with  the  highest  weight,  matching  it  with  the  system  with  the 
lowest 

m 

!>/,+ 

j- 1 

and  then  scheduling  it  with  the  crew  with  the  least  number  of  jobs  currently  in  queue.  This  is 
continued  until  all  missions  have  been  scheduled.  The  dispatching  rule  does  not  account  for  what 
subsystems  are  required  by  the  mission  in  making  the  system  assignments.  It  just  encapsulates  all 
processing  times  and  ready  times  into  a  sum  regardless  of  whether  they  will  be  required  by  the 
mission  in  the  assignment.  Table  8  presents  the  dispatching  rule’s  assignment  and  schedule  for 
our  numerical  example.  The  TWCT  objective  function  for  the  dispatching  rule  on  our  numerical 
example  is  1 1,309  which  is  41.5%  higher  than  the  objective  given  by  heuristic  solution  approach. 
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Table  8.  Dispatching  Rule  for  Numerical  Example 


Mission 

Weight  (w >i) 

m  1  m 

i'.+Jt'. 

>i  V  >=' 

System 

Crew 

ASC 

10 

2.41 

19 

1 

ASC 

10 

4 

9 

2 

ASC 

10 

6.74 

10 

3 

ASC 

10 

9.90 

21 

4 

ASN 

9 

11 

1 

5 

FSL 

8 

11 

16 

6 

FSL 

8 

12.10 

8 

1 

FSL 

8 

12.80 

5 

2 

FSL 

8 

13 

3 

3 

FSL 

8 

13.93 

13 

4 

FSL 

8 

14.86 

4 

5 

FSL 

8 

18 

11 

6 

ADC 

6 

18 

12 

1 

ADC 

6 

19.30 

23 

2 

ADC 

6 

19.97 

24 

3 

ASY 

5 

20 

14 

4 

ASY 

5 

20.97 

6 

5 

ASY 

5 

28.41 

17 

6 

ASY 

5 

28.86 

18 

1 

ASY 

5 

56.16 

15 

2 

DSP 

3 

58 

2 

3 

DSP 

3 

95.86 

7 

4 

TNG 

1 

97.41 

20 

5 

TNG 

1 

102.20 

22 

6 

24 


7.  Experimental  Design 


Realistic  problem  instances  of  the  multi-state  selective  maintenance  problem  were 
randomly  generated  and  both  the  heuristic/optimization  and  dispatching  rule  approaches  were 
tested  for  their  performance  in  terms  of  solution  quality  (as  measured  by  the  TWCT  objective) 
and  computation  time. 


7.1  Number  of  Systems,  Missions,  and  Crews 

We  will  be  using  the  F-16  as  our  motivating  example  as  outlined  above  in  section  4. 
Numerical  examples  will  be  evaluated  at  the  squadron  level  (24  planes).  All  scenarios  will  be 
assumed  to  have  6  identical  crews  available  for  maintenance  at  any  given  time.  In  each  scenario 
the  number  of  systems  is  considered  to  be  equal  to  the  number  of  missions  although  in  actuality 
the  number  of  systems  may  exceed  the  number  of  missions,  but  never  fall  below.  Thus,  there  will 
be  24  missions;  each,  one  of  the  eight  different  missions  types  outlined  above.  The  number  of 
each  mission  type  will  be  distributed  according  to  Table  9. 

Table  9.  Mission  Type  Distribution 


Mission  Type 

#  of  Missions 

FSL 

~DU[1, 24] 

ADC 

~DU[1, 24  -  #FSL] 

ASC 

~DU[1, 24  -  (#FSL  +  #ADC)] 

ASY 

~DU[1, 24  -  (#FSL  +  #ADC  +  #ASC)] 

ASN 

~DU[1, 24  -  (#FSL  +  #ADC  +  #ASC  +  #ASY)] 

DSP 

~DU[1, 24  -  (#FSL  +  #ADC  +  #ASC  +  #ASY  +  #ASN)] 

TNG 

~DU[1, 24  -  (#FSL  +  #ADC  +  #ASC  +  #ASY  +  #ASN  +  #DSP)] 

TST 

24  -  (#FSL  +  #ADC  +  #ASC  +  #ASY  +  #ASN  +  #DSP  +  #TNG) 

By  using  the  discrete  uniform  (DU)  distribution  we  ensure  that  the  number  of  missions  takes  on 
whole  number  values.  By  ordering  the  distribution  in  this  way  we  are,  in  effect,  weighting  the 
first  missions  higher  because  there  is  a  higher  likelihood  that  they  will  have  a  greater  number  of 
missions.  Whereas  the  latter  missions  (e.g.  TNG  and  TST)  will  likely  have  a  very  small  number 
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of  missions  or  none  at  all  because  all  of  the  24  missions  will  already  be  allocated  by  the  time 
these  missions  receive  their  assignments.  We  believe  this  ordering  to  be  consistent  with  the  types 
of  missions  required  by  the  Air  Force,  with  the  majority  of  missions  requiring  the  Full  System 
List  (FSL)  and  tapering  down  to  only  a  few  Training  (TNG)  and  Testing  (TST)  missions. 

7.2  State  of  the  System 

For  each  subsystem  for  a  given  system,  whether  the  subsystem  was  functioning  or  failed 
upon  returning  from  a  mission  is  considered  to  be  a  Bernoulli  trial  with  probability  of 
functioning  p  =  0.95.  Therefore,  each  subsystem  for  a  given  system  has  a  5%  chance  of  being 
failed.  Assuming  a  series  system  for  an  aircraft  with  41  subsystems,  this  would  mean  that  the  F- 
16  only  has  a  0.9541  =  0.1221  total  reliability  which  is  much  lower  than  actuality.  A  true  F-16  has 
an  estimated  reliability  of  approximately  0.90.  This  would  mean  that  each  individual  subsystem 
would  actually  have  a  0.90l/41  =  0.9974  reliability  given  a  series  system.  However,  for  the 
purpose  of  our  experiment,  a  0.9974  subsystem  reliability  is  too  high  and  does  not  generate 
enough  failures  to  test  our  solution  approaches,  therefore  0.95  was  used. 

If  a  subsystem  was  determined  to  be  failed  by  the  Bernoulli  trial,  the  number  of  hours 
required  to  bring  the  subsystem  to  a  fully  functioning  state  was  determined  by  a  lognormal 
distribution.  It  was  desired  to  have  a  distribution  whose  mean  was  4  hours  but  still  had  a  5% 
chance  of  maintenance  times  that  exceeded  24  hours.  This  was  designed  to  capture  the  fact  that 
the  average  repair  time  for  an  F-16  requires  four  hours.  Some  repairs  may  only  take  an  hour  as  a 
failed  component  may  be  easily  isolated  and  replaced.  However,  there  is  a  chance  that  an  engine 
failure  may  occur  or  a  problem  may  take  a  very  long  time  to  troubleshoot,  therefore  we  allow  for 
the  small  possibility  of  a  very  long  maintenance  time.  The  lognormal  distribution  was  chosen  for 
it  skewed  shaped  and  that  it  only  yields  positive  values.  Lognormal  distribution  parameters  with 
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an  exact  mean  of  4  and  a  Pr(X  <  24)  =  0.05  could  not  be  obtained  because  imaginary  roots  were 
obtained  when  solving  the  systems  of  equations.  The  closest  real  distribution  values  that  can  be 
obtained  are  p.  =  0.29034  o  =  1.7608  which  gives  a  mean  of  approximately  6.3.  Thus,  the  time  to 
bring  a  failed  subsystem  to  a  functioning  state  is  distributed  LN(0.29034,  3.1004).  Because  the 
lognormal  distribution  is  continuous,  and  we  desire  to  have  integer  values,  the  ceiling  of  all 
values  of  the  random  distribution  is  taken. 

7.3  Ready  Times 

Ready  times  are  assumed  to  follow  a  custom  made  discrete  distribution  as  follows 

0  hours  p  =  0.1 

2  hours  p  =  0.3 

iJ  4  hours  p  =  0.3 

96  hours  p  =  0.3 

It  is  assumed  that  if  a  spare  parts  delay  is  incurred  it  will  fall  under  one  of  the  four  scenarios 
outlined  above.  There  is  a  10%  chance  that  a  spare  part  will  not  be  needed  or  it  is  immediately 
available  on  hand  to  maintenance  personnel.  The  2-hour  spare  part  delay  scenario  is  designed  to 
model  a  part  stored  in  a  warehouse  on  the  base  that  must  be  located  and  transported  to  the 
maintenance  hangar.  The  2-hour  scenario  is  assumed  to  occur  30%  of  the  time.  The  4-hour 
scenario  models  the  cannibalization  of  a  part  from  another  plane.  It  is  assumed  to  take 
approximately  4  hours  for  the  part  to  be  extracted  from  the  cannibalized  aircraft  and  installed  in 
the  aircraft  being  maintained.  This  scenario  is  also  considered  to  occur  30%  of  the  time.  The 
final  scenario,  also  with  a  30%  probability,  is  if  the  spare  part  need  is  not  available  on-site  and 
must  be  shipped  from  a  depot  in  a  foreign  location.  It  is  believed  to  take  approximately  96  hours 
(4  days)  to  receive  a  part  from  the  depot,  and  thus  delay  maintenance  actions  for  that  amount  of 
time. 
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7.4  Weights 

The  weight  of  each  mission  was  distributed  DU[1,  10].  Thus  each  mission  was  assigned  a 
whole  number  weight  value  between  1  and  10,  inclusively.  Different  missions  were  allowed  to 
have  the  same  weight. 
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8.  Results 


There  were  233  replications  of  the  above  experimental  design  generated  and  solved  using 
both  the  heuristic/optimization  and  the  dispatching  rule  solution  approaches.  The  heuristic  coded 
in  VBA  ran  in  a  fraction  of  a  second,  however  the  optimization  of  the  scheduling  problem  using 
CPLEX  took  an  average  of  7.5  minutes  with  the  longest  solution  times  taking  over  70  minutes. 
The  dispatching  rule  runs  in  a  fraction  of  a  second  and  is  able  to  obtain  solutions  that  are  on 
average  a  mere  0.33%  above  the  heuristic/optimization  approach  and  even  beat  the 
heuristic/optimization  approach  in  86  out  of  the  233  experiments. 

It  is  impossible  to  determine  how  the  heuristic/optimization  and  dispatching  rule 
approaches  compare  to  the  optimal  solution  because  total  enumeration  is  infeasible  for  practical 
size  problem  instances.  It  is  reasonable  to  conclude  that  both  solutions  approaches  are  providing 
quality  solutions  due  to  the  similarity  of  the  results. 
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Appendix  A 


File  name:  MSSMx_Solution_Code.xls 
1 .  Setting  up  the  spreadsheet  for  use 

When  you  open  the  file,  you  may  receive  a  message  similar  to  the  one  shown  in  Figure 


Figure  A.  1  Macro  Notification 


The  Visual  Basic  code  used  to  evaluate  the  model  is  written  within  macros.  Therefore,  you 
should  click  on  “Enable  Macros.”  This  will  open  the  “missions”  worksheet  shown  in  Figure  A.2. 
2.  Inputting  Data 

This  “missions”  worksheet  (Figure  A.2)  displays  the  subsystems  required  by  each 
mission.  Next,  enter  the  number  of  crews  available  to  perform  maintenance  at  any  given  time. 
Then  enter  the  number  of  missions  of  each  type  needed  to  be  performed  in  the  row  of  the 
corresponding  three-letter  mission  abbreviation  in  the  field  “#  Missions”.  Likewise,  enter  the 
weight  or  importance  of  the  mission  in  the  “Weights”  field.  The  weight  should  be  a  whole 
number. 
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Figure  A.2  Missions  Worksheet 


To  enter  number  of  systems  available  and  the  status  of  each  system,  click  on  the  a_ij  worksheet 
tab  to  display  the  a_ij  worksheet  shown  in  Figure  A.3.  Enter  the  number  of  systems  available  in 
the  “Number  of  Systems  Field.”  Next,  enter  the  number  of  maintenance  hours  required  to  bring 
the  subsystem  of  a  given  system  to  a  fully  functioning  state.  This  should  be  a  whole  number 
value.  Match  the  column  of  the  corresponding  system  number  with  the  row  of  the  subsystem 
requiring  maintenance. 
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Figure  A.3  a_ij  Worksheet 


Next,  click  on  the  rho_ij  tab  at  the  bottom  of  the  screen  to  display  the  rhoij  worksheet  show  in 
Figure  A.4.  This  worksheet  will  allow  the  user  to  input  the  ready  times  of  the  subsystems  for 
each  mission.  Ready  times  should  be  whole  number  values  and  entered  exactly  as  the 
maintenance  hours  described  above  for  worksheet  a  ij. 
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Figure  A.4  rho_ij  Worksheet 


Once  all  of  the  system  input  parameters  have  been  entered,  return  to  the  missions  worksheet  by 
clicking  on  the  “missions”  tab  at  the  bottom  of  the  screen.  This  will  return  the  user  to  the  screen 
shown  in  Figure  A.2.  Click  the  button  labeled  “Heuristic”  in  the  upper  left-hand  comer  of  the 
screen  to  run  both  the  heuristic  and  dispatching  rule.  The  assignment  made  by  the  heuristic  will 
be  written  to  the  “assign”  worksheet  as  show  in  Figure  A.5.  This  format  is  the  exact  input  format 
required  by  AMPL  to  use  CPLEX  to  solve  for  the  schedule.  This  text  may  be  pasted  into  a  text 
file  to  be  read  in  by  AMPL.  AMPL  code  is  also  contained  in  this  appendix. 
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Figure  A.5  Assign  Worksheet 

The  results  of  the  dispatching  rule  are  written  to  “dispatch”  worksheet  show  in  Figure  A.6.  The 
dispatch  worksheet  shows  the  in  descending  order,  going  down  the  rows,  the  order  in  which  the 
mission/system  assignments  are  dispatched.  The  “k”  column  refers  to  the  mission  number  and 
the  “i”  column  refers  to  the  system  number.  Mission  and  system  numbers  on  the  same  row  are 
assigned  to  each  other.  The  “w”  column  displays  the  weight  of  the  mission  and  the  ratio  column 

m  m 

shows  the  ^atJ  +  calculation.  The  objective  function  of  the  dispatching  rule  is  displayed 

M  V  j= i 


on  the  upper  right-hand  side  with  a  “1”  next  to  it  indicating  the  first  solution  in  the  experiment. 


36 


Figure  A.6  Dispatch  Worksheet 
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param  gamma; 
param  k; 


AMPL  Scheduling  Formulation 


set  crew  =  1  ..  gamma; 
set  mission  =  1  ..  k; 

param  w  {mission}; 
param  p  {mission}; 
param  r  {mission}; 

param  cmax  =  sum  {j  in  mission}  (r[j]  +  p[j]); 

set  time  =  0  ..  cmax-1; 

var  s  {mission}  >=0  integer; 
var  x  {crew,  mission,  time}  binary; 

minimize  TWCT:  sum  {j  in  mission}  w[j]  *  sum{i  in  crew,  t  in  time}  (t  +  p[j])  *  x[ij,t]; 
subject  to  one  {j  in  mission}:  sum  {i  in  crew,  t  in  time}  x[ij,t]  =  1; 

subject  to  two  {i  in  crew,  t  in  time} :  sum  {j  in  mission,  tt  in  max(t-p[j],0)  ..  t-1 }  x[i J,tt]  <=  1 ; 
subject  to  three  {j  in  mission}:  sum{i  in  crew,  t  in  time}  t  *  x[ij,t]=  s[j]; 
subject  to  four  {j  in  mission}:  s[j]  >=  r[j]; 
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